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Abstract. The Casimir effect in quantum electrodynamics (QED) is perhaps the best-known 
example of fluctuation-induced long-ranged force acting on objects (conducting plates) immersed 
in a fluctuating medium (quantum electromagnetic fleld in vacuum) . A similar effect emerges in 
statistical physics, where the force acting, e.g., on colloidal particles immersed in a binary liquid 
mixture is affected by the classical thermal fluctuations occurring in the surrounding medium. 
The resulting Casimir-like force acquires universal features upon approaching a critical point 
of the medium and becomes long-ranged at criticality. In turn, this universality allows one 
to investigate theoretically the temperature dependence of the force via representative models 
and to stringently test the corresponding predictions in experiments. In contrast to QED, the 
Casimir force resulting from critical fluctuations can be easily tuned with respect to strength 
and sign by surface treatments and temperature control. We present some recent advances 
in the theoretical study of the universal properties of the critical Casimir force arising in thin 
films. The corresponding predictions compare very well with the experimental results obtained 
for wetting layers of various fluids. We discuss how the Casimir force between a colloidal particle 
and a planar wall immersed in a binary liquid mixture has been measured with femto-Newton 
accuracy, comparing these experimental results with the corresponding theoretical predictions. 



1. Introduction 

This contribution focuses on the Casimir effect which arises when one confines fiuctuations 
of a different nature compared to those originally considered by Casimir in his pioneering 
work [Ij. Indeed we shall be concerned here with the fluctuations of thermal origin which 
occur close to a second-order phase transition point {critical point). The resulting Casimir effect 
is apparently less widely known than the corresponding one within quantum electrodynamics 
(QED). Nonetheless, the resulting femto-Newton forces at the sub-micrometer scale can be put 
to work in soft matter systems and their high degree of tunability might even be exploited for 
concrete applications in colloidal suspensions. However, before focussing on the Casimir effect 
due to critical fluctuations we first discuss in subsec. 12.11 the effect within QED, following in 
spirit the original derivation by Casimir [T]. In subsec. [2^2] we then explain how this effect arises 
in statistical physics upon approaching a critical point, highlighting the analogies and differences 
between the two. In sec. Owe present an overview of the available theoretical predictions which 
are relevant for the qualitative and quantitative interpretation of the experimental results. A 
short summary, with perspectives and applications is presented in sec. HI 



2. The Casimir effect 

2.1. The effect in Quantum Electrodynamics 

The Casimir effect is named after Hendrik Casimir who discovered, 60 years ago, that - quite 
surprisingly - two parallel, perfectly conducting and uncharged metallic plates in vacuum attract 
each other due to the quantum fluctuations of the electromagnetic fields, even at zero temperature 
r = [Ij. Indeed the plates (assumed to be perfectly conducting) effectively impose boundary 
conditions (BCs) on the electromagnetic fields so that, more specifically, Ey, = where 
E|| and are the components of the electric E and magnetic B fields which are parallel and 
transverse to the surface of the plates. As a result of these boundary conditions, the fluctuation 
modes of the fields in the space within the two plates can only have a specific set of L-dependent 
allowed wave-vectors, where L is the separation between the parallel plates. For example, in the 
case of Dirichlet BCs, the component k± of the wave- vector k perpendicular to the plates assumes 
only the quantized values k±_ = TmjL with n = 1,2,.... Roughly speaking, the "unbalance" 
between the pressure exerted by the allowed modes within the plates and the one exerted by 
the modes outside them is at the origin of the Casimir effect. This statement can be made 
quantitative by calculating the size-dependent energy of the electromagnetic fields allowed 
in the portion of the vacuum within the plates, which is given, according to QED, by 



£{L) = ^ ^?lc|kmodes| , 
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where kmodes are the corresponding wave-vectors and c is the speed of light. For perfectly 
conducting plates of large transverse area S, the expansion of £{L) in decreasing powers of L 
takes the form (in three spatial dimensions) 
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where £'buik oc SL^ is the energy associated to the electromagnetic field of the vacuum in 
the absence of the plates within the volume of space V = SL enclosed by them; £^1^^^^^ = 

^surf + ^surf ^ °^ energies £suvi associated to the introduction of each single 

plate, separately, in the vacuum; The next term in the expansion i.e.. 
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is proportional to SL~^ and represents the interaction energy between the two plates which is 
due to their simultaneous presence in space. This is the Casimir term we are interested in and 
which is responsible for the celebrated Casimir force. Further terms in the expansion are of order 
0((kL)~^) compared to £cas, where k is a material-dependent parameter which describes the 
deviations of the plates from the perfectly conducting behaviour due to the finite conductivity of 
the metal. In the expansion ([2]) the first term £^buik is a property of the vacuum, the second one 

(^surf^^) depends on the material properties of the plates, whereas the third one, £'cas in eq. 1^, 
is universal in that it does not depend on the specific material which the plates are made of, but 
only on geometrical properties {S and L) and on fundamental constants (he). The relevance of 
higher-order terms in the expansion ^ depends on the material parameter k. 

The Casimir term, the third of the the expansion of the vacuum energy £{L) in decreasing 
powers of L, actually gives rise to a measurable effect: a small displacement 6L of one of the 
two plates results in a change 5£{L) in the energy of the fields within the plates and therefore 
in a force Fi^i{L) = —5£/5L (pushing on the inner surface of each single plate). Focussing on 



one of the two plates, the contribution to the force Fin(i) due to the change of £'buik in £{L), 
i.e., -Fbuik = —^£hnik/SL, is actually canceled out by the force Fout = F\a{L = oo) = Fbuik 
acting from the other side of the plate and due to the fluctuation modes in the corresponding 
half-space outside the plate. £'surf does not change upon changing L so that the total force 
F[L) = Fin(-L) — Fout acting on each single plate is due only to the change in the Casimir 
energy: 

The universal behaviour encoded in this equation can only be detected for a separation L 
between the plates which is much larger than the length scale being material-specific 

properties relevant at smaller distances. In particular, k is related to the typical wave-vector 
at which the plates are no longer effective in imposing the boundary conditions on the fields. 
As a rough estimate, this scale is set by the plasma frequency ujp of the metal (more detailed 
calculations support this argument p]), i.e., ck ~ cOp, where cjp ~ 3 • lO^'^^Hz for copper, yielding 
~ 0.3/im. Attempts to verify experimentally eq. (jl]) started already in 1958 [3J but the 
first sound experimental confirmation came only 40 years later within the range of separations 
0.5/xm ^ L < 6/im [2, 4J. Apart from the relevance of this effect for possible technological 
applications in micro- and nano-electromechanical systems (MEMS and NEMS, respectively), 
in which the associated force is responsible for stiction, its precise measurement could be useful 
in order to test of the validity of some fundamental laws down to the nano- meter scale, see, e.g., 
ref. [5]. 



2.2. The effeet in Statistieal Physics 

Thirty years after the seminal paper by Casimir, M. E. Fisher and P. -G. de Gennes published 
a note "On the phenomena at the walls in a critical binary mixture" [6] in which it was shown 
that Casimir- like effects (i.e., fluetuation-induced forces) arise also in statistical physics when a 
medium in which fluctuations of a certain nature take place is spatially confined |7J . In order to 
understand how these effects arise and what their relation is with the Casimir effect described in 
subsec. 12.11 we briefly recall here some basics facts about binary mixtures, their critical points 
and the effect of confinement. The schematic phase diagram of a liquid binary mixture in the 
bulk is depicted in fig. [H For a fixed pressure, the relevant thermodynamic variables are the 
temperature T and the mass fraction ca of one of the two components in the mixture. For a 
certain concentration ca, the two components of the mixture are mixed at low temperatures and 
the resulting liquid solution looks homogeneous in a test tube (schematically represented on the 
left of the phase diagram). Upon increasing the temperature, however, the liquid demixes into an 
A- and a B-rich solution and becomes inhomogeneous in the test tube (represented on the right), 
where the two solutions are typically separated by an interface. The transition between the mixed 
and the demixed phases occurs at the solid first-order transition line. At the lowest point (CP) of 
this line the transition becomes critical, i.e., second-order. This point is referred to as the (lower) 
critical point. The mixed and demixed phases can be distinguished by monitoring the so-called 
order parameter of the transition, which can be identified with the deviation 5ca(x) = ca{x) — ca 
of the local concentration ca(x) of, say, component A from its average value in the mixture. 
Indeed, the thermal average ((5c^(x)) of the order parameter is homogeneous in the mixed 
phase, whereas it assumes two different values in the A-rich and B-rich solution which coexist 
in the demixed phase. The order parameter is characterized by a spatial correlation length ^, so 
that, within each single phase, {6ca{'x)Sca{x')) — (5ca(x))(5ca(x')) oc exp{— |x — x'|/^}, where ^ 
depends on T and ca- Upon crossing the first-order transition line, the correlation length stays 
finite, whereas it diverges upon approaching the critical point. The physical behaviour of the 
system at the scale set by is actually determined by the fluctuations of the order parameter, 
which becomes the relevant physical quantity for the description of the system. In terms of it. 




Figure 1. Schematic phase diagram of a binary hquid mixture with a lower critical point (CP). 
On the left and on the right of the diagram, the schematic side views of a test tube filled with 
the binary liquid mixture respectively in the mixed and in the demixed phase are shown. 



one can express the effective free energy T of the system and determine the behaviour of the 
thermodynamic properties upon approaching the transition. 

In order to understand the consequences of confining such a binary liquid mixture, we have 
first to consider the result of inserting a single plate (made up of, say, glass) into the mixture. In 
this case, due to expected quantitative differences in the interactions between the molecules of 
the plate substrate and those of the two components of the mixture, the plate shows preferential 
adsorption for one of the two components, say, A. Accordingly, the presence of the plate in 
the mixture induces a local increase of the average order parameter ((5ca(x)) close to it. Upon 
increasing the correlation length ^ (i.e., getting closer to the critical point), the plate effectively 
imposes boundary conditions on the order parameter [8l|9]. 



Table 1. Analogies and differences between the Casimir effect in (zero-temperature) QED and 
in Statistical Physics (see also footnote [I]) . 



fluctuating quantity: 
excitation: 

range of fluct.: 



QED 
E, B 

Quantum 
he (T = 0) 

oo 



Stat. Phys 
order parameter 

Thermal(classical) 
k^T {h = 0) 

finite: ^ 

^ / oo close to CP 



Confinement 



4 

long-ranged force 



range: ^ 
long-ranged at CP 



Now that we have recalled how to describe the binary mixture in terms of a fluctuating order 
parameter (generically refereed to as (j) in what follows) and the effect of plates in terms of 
boundary conditions, it is easy to understand in which sense a Casimir-like interaction arises in 
these instances. Indeed (see tabled]), the Casimir effect in QED emerges because there are some 
quantities (the fields E and B) which "fluctuate" in space and time due to the quantum nature 
of the "medium" (actually the vacuum) in which their fluctuations take place. Accordingly, the 
relevant scale of the phenomenon is set by he at zero temperature. The two-point correlation 
function of the fluctuations of these quantities are characterized by an algebraic decay in space, 
i.e., the associated range is infinite (due to the vanishing mass = of the photon). When 
the spectrum of these fluctuations gets modified by external bodies (metallic plates) imposing 
boundary conditions on the fluctuating fields, the associated energy £ turns into a function of 
the position of these boundaries and as a result an effective long-ranged Casimir force arises on 
them. In statistical physics, close to a critical point, the relevant fluctuating quantity is the order 
parameter <j) of the phase transition (</>(x) = (5c^(x) in the case of the binary mixture) and its 
fluctuations are of thermal nature and due to the coupling to the thermal bath at temperature 
T. Accordingly, the relevant scale of the phenomenon is k-^T and the fluctuations are of classical 
nature (i.e., h plays no role). The range of these fluctuations is given by ^ and therefore is finite 
but tunable upon approaching the critical point. As in the case of QED, external bodies impose 
boundary conditions on the order parameter (/) and therefore they affect the spectrum of its 
allowed fluctuations and the associated effective free energy T depends on the positions of these 
bodies. As a result, a force is expected to act on them, with a range set by the correlation length 
^. In table [1] we summarize analogies and differences between the Casimir effect occurring in 
QED and in statistical physical. Even though we have been referring to the case of a binary 
liquid mixture, the line of argument presented above extends to critical phenomena in general. 
Indeed, consider a medium close to its critical point and confined between two parallel plates 
at a distance L (film geometry). The medium might be, e.g., pure ^He or a ^He-^He mixture 
close to the superfluid transition, a binary liquid mixture, a Bose gas close to the Bose-Einstein 
condensation, liquid crystals etc. If the correlation length of the fluctuations of the order 
parameter and the thickness L of the film are much larger than the microscopic length scale 
^micr (set, Say, by the molecular scale of the system), the free energy T of the confined medium 
can be decomposed in decreasing powers of L for a fixed value of L and large transverse area 
5", in analogy to eq. (f2]) [TT] 

.^(T, L) = ^buik + -^Lf^ + S^e,! (L/0 + corr. (5) 

We have assumed here that the only relevant thermodynamic variable is the temperature 
r, which determines the correlation length ^ ~ '?o|(2^ ~ T^c)lT^c\~^ where is a system- 
specific quantity and v a universal exponent, in the sense specified below. Possible additional 
thermodynamic variables, such as the concentration ca of the binary mixture, are assumed to be 
tuned to their critical values, i.e., ca = ca,c with reference to fig. E In the expansion ([5]) the first 
term ^buik is the free energy of the bulk medium in a volume S x L, i.e., ^buik oc SL, the second 
one i^surf^^) represents the sums of the free energy costs T^^l^^ oc S for the separate introduction 

of the two walls in the system, i.e., J^^^^^f^^ = J^sulf + -^iurf "^^^ third term represents, as in 
eq. the interaction (free) energy between the two walls. Compared to the analogous term 
in the case of QED we note, as expected, that k^T replaces he in providing the scale of the 
phenomenon d and that the interaction between the two walls is no longer long-ranged but its 

^ In some cases, the correspondence illustrated in table[l]goes beyond the mere analogy. Indeed, for the simplest 
geometrical settings, such as film geometry, there is a mapping between the Casimir effect in QED in d spatial 
dimensions and the effect at the critical point of the so-called Gaussian model in d + 1 dimensions. 
^ Due to the different engineering dimensions of [he] = energy x length and [ksT] — energy, and taking into 



range is set by ^ via the dependence on L/^ of the scahng function 0||(p). This dependence can 
be understood as follows: the perturbation that the wall induces on the order parameter extends 
within a distance ~ ^ far from it, so that the two walls can interact with each other only at 
separation L < ^ whereas the interaction energy vanishes for L S> ^, so that (x oo) = 0. 

The scaling function Gy is universal (as the analogous term £'cas dS]) in eq. ([2])) in the sense 
that it depends only on some gross features of: 

(a) the system in the bulk, such as the range and symmetries of the interaction and the kind 
of order parameter which describes the phase transitions, 

(b) the surfaces which provide the confinement, such as the possible symmetries of the bulk 
system which they break, e.g., by favoring certain values of the order parameter at the 
boundaries. These preferences eventually translate into effective boundary conditions for 
the order parameter [8l|9]. In addition. By depends on the shape of the boundaries, as we 
shall see below. 

The general features at points (a) and (b) define the so-called bulk and surface universality classes 
of the confined system, respectively. This universality is typical of critical phenomena and allows 
one to investigate universal properties such as i' and Gy by means of suitable representative 
models belonging to the same bulk and surface universality classes of the actual system one is 
interested in and which lend themselves for a simpler theoretical analysis. Analogously to what 
happens in QED, a small change SL of the distance between the confining surfaces leads to a 
change 5 J- in the free energy and therefore to a force —6J-/5L acting on the displaced wall. The 
L-independent contribution to this force due to .^i^uik in eq. ([5]) is counterbalanced by the same 
contribution acting from outside the walls when they are immersed in the critical medium, so 
that the net force F{L,T) is given by (see eq. ([5])) 

E^^ = t^^^^{L/0 for L,^»W, (6) 

where ^\\{u) = — 2G|[(n) + uGj|(ti) and it shares the properties of universality with Q(u). The 
universal scaling function can be conveniently determined from the analysis of suitable 
representative models, either via field-theoretical methods or Monte Carlo simulations. In turn, 
the resulting predictions can be very stringently compared with the experimental determination 

of 

2. 3. The force of quantum and critical fluctuations 

In the preceding two sections we have briefly reviewed the origin of the Casimir effects due 
to quantum and critical fluctuations, highlighting the analogy and differences between them 
and the corresponding scaling and universal properties. In view of possible experimental 
investigations, however, it is important to estimate the expected magnitude of these effects. 
Assuming (quite unrealistically, indeed) that one is able to realize a system of two parallel plates 
of area S = 1 cui^ at a distance L = 1 //m, the zero-temperature Casimir force resulting from 
quantum fluctuations is roughly estimated as ~ x hc/L'^ ~ 6 • 10~^N, whereas the force due 
to critical fluctuations possibly occurring at room temperature (T ~ 300 K) can be estimated as 
~ 5 X k-QT/Lp' ~ 4 • 10~^ N. Even though these appear to be quite small forces, their magnitudes 
are actually comparable with the weight ~ 2 • 10~^ N of a water droplet of half a millimeter in 

account the expected proportionality to 5", the two effects are characterized by different powers in their dependence 
on L. 

^ The scahng function By actually depends also on whether a certain correlation length ^ is realized above or 
below the bulk critical temperature Tc- This is made clear later on by expressing the scaling function as a function 
of an appropriate scaling variable x, c.f., eq. instead of u = L/^, as done here. 




Figure 2. Formation of a complete wetting film of a fluid. Close to the liquid-gas phase transition 
(a), a vapour (b) condenses in the presence of a suitably chosen (solid) substrate, forming a liquid 
film (c), the thickness L of which diverges upon approaching the bulk condensation pressure Pq 
(d). 



diameter. However, the realization of such a geometrical setting is actually extremely difficult 
due to the problem of maintaining the alignment between the two plates within the required 
accuracy. This problem is particularly difficult when trying to detect the Casimir effect in 
quantum electrodynamics and indeed this has been achieved in the parallel-plate geometry only 
quite recently [4j . For the detection of the critical Casimir effect the problem of alignment can 
be solved by taking advantage of physical phenomena, such as wetting, which naturally lead to 
the formation of confined films of fluids of a certain controllable well-defined thickness L, which 
can be varied by acting on the thermodynamic parameters of the system, such as pressure and 
temperature [l2l [13] . 

2.4- Wetting films and critical endpoints 

In fig. we depict schematically the phase diagram of a substance, as a function of the 
pressure P and temperature T, showing the generic solid (sol), liquid (liq), and gas phases. In 
order to obtain a liquid wetting film of this substance, the temperature Tq and pressure P have 
to be chosen in such a way that the bulk system is in the vapour phase (fig. [^^b)) but close to 
the condensation transition, occurring at P = Po{Tq). If the vapour is taken in contact with a 
suitably chosen (solid) substrate (hatched grey area in fig. [2l^c)), also at temperature Tq, a fiuid 
film will condense on it as a consequence of the interaction of the substance with the substrate 
which makes the formation of a liquid layer thermodynamically favourable. The thickness L 
of the resulting film can be controlled acting on the undersaturation 5P = P — Po(2o) and is 
determined by van-der-Waals and dispersion forces. Depending on the choice of the substrate, L 
diverges smoothly upon approaching the bulk condensation pressure, i.e., for 6P 0. This case 
is referred to as complete wetting and is characterized by the formation of a liquid film which can 
be made macroscopically thick, so that L ^ imicr- The phenomenon of wetting leads naturally 
to the formation of a liquid film of constant thickness L, in which the liquid is confined between 
the solid substrate and the liquid-vapour interface. The latter is macroscopically well defined 
sufficiently far from the liquid-vapour critical point. Even though the problem of the alignment 
between the confining surfaces is seemingly solved by wetting films, one has to be aware of 
the fact that the liquid- vapour interface actually fluctuates {capillary fluctuations) around its 
average position. 

A possible indirect evidence of critical Casimir forces can be obtained by monitoring the 
thickness L of the wetting layer upon approaching a critical point which, in the bulk, occurs 



within the hquid phase, is close to the Uquid-vapour first-order transition but is far enough from 
the hquid- vapour critical point. This is the case for critical endpoints (cep), which are located 
where the liquid-vapour transition line terminates a line of critical points (in the liquid phase, 
see, c.f., fig. El^a)). Indeed, upon approaching a critical endpoint, the critical fluctuations of 
the associated order parameter, confined within the liquid wetting film of thickness L, give rise 
to the critical Casimir effect. The associated force adds up to the previously acting (van-der- 
Waals) forces in determining the equilibrium distance L of the liquid-vapour interface from the 
substrate and consequently it affects the thickness of the wetting film. The dependence of L 
on the thermodynamic control parameters such as pressure and temperature is the basis for 
the determination of the Casimir force within this approach, which was proposed and discussed 
in detail in ref. p[^, providing the theoretical motivation for the experimental investigations 
summarized in subsec. 13.11 

3. The critical Casimir effect at work 

In this section we summarize the theoretical predictions for the scaling function of the critical 
Casimir force in thin films and we compare them with the available experimental results based on 
wetting films (subsec. l3. ID . Then we discuss a more direct measurement of this force (subsec. l3.2p . 

3. 1 . Wetting films 

In subsec. 12.41 we argued that it is possible to exploit critical endpoints in order to infer the 
critical Casimir force from the thickness of a complete wetting film. One of the most studied 
fluids with such a critical endpoint is pure ^He, the phase diagram of which is sketched in fig.[3]^a). 
Within the fluid phase ^He undergoes a second-order phase transition between a superfluid (also 
refereed to as Hell) and a normal (Hel) behaviour upon increasing the temperature across the 
so-called A-line. The order parameter (j) of this phase transition is physically provided by the 
wave-function of the superfluid. The A-line terminates at the critical endpoint (cep) on the fluid- 
vapour first-order transition line (which ends on the right at the liquid-vapour critical point) 
located at low pressure, as indicated by the closer view of fig. [3]^a). As previously described, in 
order to detect critical Casimir forces one prepares the system in the vapour phase corresponding 
to the point P in the closer view of fig. El^a) and exposes this vapour to a suitable substrate 
(e.g., Cu) in order to obtain a complete wetting film of thickness L. Then one changes the 
thermodynamic parameters in such a way to follow the thermodynamic path indicated as 7 
in fig. E^a), along which the film thickness L would not change in the absence of the critical 
endpoint. The actual changes are therefore due to the action of the critical Casimir force on the 
liquid-vapour interface and from these changes it is possible to determine experimentally the 
scaling function of the force (see eq. ([6])) [13J. In order to predict theoretically the scaling 
function one has to single out a theoretical model which belongs to the same bulk and surface 
universality class as the confined ^He. It is well known that the universal aspects of the bulk 
critical behaviour of ^He are properly captured, upon approaching the A-line, by the lattice XY 
model [H], which lends itself, e.g., for Monte Carlo simulations. Accordingly, the behaviour of 
a film of ^He close to the A-line is captured by the XY model with film geometry, where the 
boundary conditions on the lattice are chosen to be in the same surface universality class (i.e., to 
lead to the same effective boundary conditions for the order parameter) as the actual confining 
boundaries, represented by the substrate-liquid and liquid-vapour interfaces. As previously 
mentioned, the order parameter <j) of the superfluid transition is related to the wave- function 
of the superfluid component and therefore it has to be continuous in space and to vanish both 
inside the substrate and within the vapour. As a result, it has to vanish at the boundaries, i.e., 
(j) satisfies Dirichlet BCs at both surfaces (DD). On the lattice model these are realized via free 
boundary conditions for the degrees of freedom. Having identified a model which belongs to the 
same bulk and surface universality class as the film of ^He, one can proceed to the determination 




Figure 3. (a) Phase diagram of pure ^He. The closer view highlights the presence of the critical 
endpoint (cep). (b) Universal scaling function •d^^ of the critical Casimir force within the three- 
dimensional XY universality class and Dirichlet-Dirichlet (DD) BCs, as a function of the scaling 
variable x = {T/Tx - l)(L/^o)^''''- The data points were obtained by Monte Carlo simulation 
of the XY model on the lattice in ref. [Ej, whereas the solid line represents the experimental 
results of ref. [19] . 



of the scaling function "i^jj . In fig. ^h) the data points refer to the scaling function ■d^^ as inferred 
from Monte Carlo simulations of the XY model in film geometry of various thicknesses and DD 
boundary conditions [15] (see also ref. [E]). The scaling function is plotted as a function of the 
proper scaling variable 

X = t{L/i,f'- (7) 

where t = (T — T\)/Tx, T is the temperature at which the lattice model is simulated and Tx is 
the associated bulk critical value (corresponding to the A-line in the actual system). L is the 
thickness of the film, the non-universal amplitude which controls the divergence of the bulk 
correlation length ^ and which can be determined independently via Monte Carlo simulation, 
and v ~ 0.66 [14J is the universal critical exponent of the correlation length. The scaling variable 
X in eq. d?]) is related to u = L/^^ in eq. Q by x = u^^'^ for t > 0. As a confirmation of the scaling 
behaviour in eq. Q, the scaling function in fig. [3]^b) does not depend on the thickness L of 
the film in which it has been numerically determined (as long as L is large enough compared 
to the lattice spacing) and the data sets corresponding to different sizes actually fall onto the 
same master curve. The predicted critical Casimir force turns out to be attractive in the whole 
range of temperatures (t?|| < 0), as expected it decays to zero for ^, i.e., for x ^ 1, whereas 
it saturates to a non- vanishing value for x ^ — oo due to the fact that at low temperatures 
long-range correlations are maintained by Goldstone modes. Some of the qualitative features 
of the present force, such as the occurrence of a deep minimum below Tx, can be understood 
within a simple mean-field approximation (see, e.g., ref. |17[ I18j). 

The very delicate experimental determination of the critical Casimir force acting on wetting 
films of ^He on a Cu-substrate has been done in ref. [19] , where the thickness L ~ 200 . . . 300 A of 
the film as a function of the distance from the critical point has been determined by capacitance 
measurements. The corresponding experimental data are reported as a solid line in fig. [3]^b) 
and show a remarkably good agreement with the corresponding theoretical predictions based 



on Monte Carlo simulations of the lattice model. In determining the value of the abscissa x 
corresponding to a certain experimental point one has to use for the system-specific value 
which has been determined independently in the experiment on the basis of the behaviour of 
the bulk correlation length. If this is properly done, then there are no free parameters which 
can be adjusted in the comparison. In this sense the test of the theoretical predictions is really 
stringent. 

On the basis of the values of the experimental parameters of ref. [K] one can estimate the 
critical Casimir pressure to be of the order of 2 Pa which, however, is still of quite difficult 
detection in those experimental conditions. In order to increase the magnitude of the critical 
Casimir force, which is proportional to the temperature T, it is convenient to look for critical 
endpoints occurring at higher temperatures compared to the one of ^He, for example by 
considering wetting films of classical binary liquid mixture. It turns out that in these cases 
the boundary conditions can be tuned in such a way to result in repulsive forces. These forces 
could find application for compensating the attractive electro dynamic Casimir force (the sign of 
which is quite difficult to be controlled) and therefore avoid stiction in micro- and nano-machines 
(MEMS and NEMS). 

The phase diagram of a binary liquid mixture is depicted in fig. [1] for a fixed value of the 
pressure P (i.e., fig. [T] is a cut of the phase diagram in the {ca,T, P)-space). Upon decreasing 
the pressure the location of the critical point CP changes in the (c^i, r)-plane and eventually 
it meets at the critical endpoint the sheet corresponding to the transition between the liquid 
and the vapour of the mixture. This critical endpoint can be exploited in order to create a 
wetting film of the binary mixture [13j, similarly to the case of ^He. The scaling function "i^n can 
therefore be inferred by monitoring the equilibrium thickness of the wetting layer as a function 
of the thermodynamic parameters. Also in this case, theoretical predictions for the scaling 
function -i^n can be obtained by studying a suitable model which belongs to the same bulk and 
surface universality class as the confined classical binary mixture. It is well-known that the bulk 
critical behaviour close to the demixing critical point in a binary mixture is properly captured 
by the lattice Ising model [H], which lends itself, e.g., for numerical studies via Monte Carlo 
simulations. The universal behaviour of the binary mixture confined in the film is therefore 
captured by the Ising model in a film geometry, in which the boundaries are chosen to be in 
the same surface universality class as the actual confining surfaces constituted by the substrate- 
liquid and liquid-vapour interfaces. In subsec. 12.21 we have already mentioned the fact that 
the boundaries generically show preferential adsorption for one of the two components of the 
mixture. As a result, the order parameter (5cyi(x) of the binary mixture (see the discussion in 
subsec. 12. 2p either increases (+) or decreases (— ) upon approaching the boundaries and this 
tendency turns into effective boundary conditions close enough to the critical point [HI E]. We 
shall refer to the case in which both boundaries preferentially adsorb the same component of 
the mixture as (++) and to the case of opposite preferences as (+— )• Additional details of 
the strength of this preferences etc. turn out to be irrelevant sufficiently close to the critical 
point. In the lattice Ising model, these boundary conditions are realized by fixing the boundary 
degrees of freedom (±1 spins) either to the same [(++)] or to opposite [(H — )] values on the two 
opposing boundaries. The resulting numerical prediction of the scaling function ??|| is reported 
in fig. [3] as a function of the scaling variable x in eq. ([7|) , where L is the thickness of the film, 
v ~ 0.63 [14J is the universal critical exponent of the correlation length, t = (T — Tc)/Tc where 
T is the temperature at which the lattice model is simulated and Tc its critical value. As 
it was the case for ^He, the value of the non-universal amplitude for this lattice model is 
measured independently via Monte Carlo simulations. The collapse of data points referring to 
lattices of different thicknesses onto the same master curve is again a confirmation of the scaling 
behaviour expected on the basis of eq. ([6]E|. The sign of the resulting Casimir force depends 

* In order to extract in this case the asymptotic scaUng behaviour from Monte Carlo data it is necessary in to 




Figure 4. Universal scaling function "i^y of the critical Casimir force within the three-dimensional 
Ising universality class and (++)/(H — ) BCs, as a function of the scaling variable x = t{L/^oy^'^ . 
The data points were obtained in ref. p5] by Monte Carlo simulation of the Ising model on the 
lattice, whereas experimental data for the case of (+— ) BCs were obtained in ref. |21] . 



on the boundary conditions and is attractive for (++) and repulsive for (H — ), corresponding, 
respectively, to < and > in the whole range of temperatures. As expected, 

these scaling functions decay to zero for x ±oo due to the fact that away from the critical 
point the correlation length becomes much smaller than the film thickness L both in the mixed 
(disordered in terms of the Ising model) and demixed (ordered) phase. In addition, the typical 
magnitude of the repulsive force for (+— ) BCs is larger than the one of the attractive force for 
(+-1-) boundary conditions. This is due to the fact that in the former case the fluctuations of the 
position of the interface between the region with positive and negative values of {6ca{x)) also 
contribute to the force acting on the confining surfaces. Also in this case some of the qualitative 
features of the force are already captured properly within a mean-field approximation [20j 

The critical Casimir effect in wetting films of classical binary mixtures has been investigated 
in refs. |2H [^5] where either ellipsometry [22] or X-ray scattering [21] techniques have been 
employed in order to determine the equilibrium thickness L of the wetting films. In particular in 
ref. [21] a binary mixture of methylcyclohexane (MC) and perfiuoromethylcyclohexane (PFMC) 
with an upper critical point at temperature Tc = 42.6°C and at a critical molar fraction of 
PFMC cpFMC,c = 0.36 has been used as the critical medium. The investigated wetting films 
were formed on a Si02/Si substrate and they had a typical thickness L ~ 100 A with the 
PFMC and MC preferring to be close to the liquid-vapour (+) and the liquid-substrate (— ) 
interface, respectively. The corresponding experimental data for the inferred scaling function of 
the critical Casimir force are reported on the right panel of fig. [H As in the comparison between 
theoretical and experimental predictions for ^He, the value of the abscissa x corresponding to 
a certain experimental point has to be calculated according to eq. ([7]) taking into account the 
experimentally determined value of ^-nd the fact that t = {T — Tc)/Tc for an upper critical 
point. Note that, again, there are no adjustable parameters for the comparison. The agreement 
with the corresponding theoretical predictions is rather good (for a more detailed discussion, see 
ref. [E]). 

All the evidence previously discussed for critical Casimir forces relies on the use of wetting 



account for corrections to scaling, see ref. [15] for details. 



films of (classical and quantum) fluids, a rather indirect way of detecting the effects of these 
forces. In addition, the reliable quantitative interpretation of the observed variations of the 
thickness of these films in terms of the action of Casimir forces requires an independent and quite 
detailed knowledge of some material-dependent parameters which determine the van-der-Waals 
forces acting within the system. In view of these potential limitations, it would be desirable to 
measure directly the critical Casimir forces, e.g., by determining the associated potentials. 

3.2. Brownian motion of colloids 

The experimental evidence summarized in the previous subsection suggests that critical Casimir 
forces become relevant at the submicrometer scale, so that, on dimensional ground and for 
critical points at room temperature Tc — 300 K one estimates the typical scale of force as 
~ k^Tc/ {0.1 fim.) — 40 fN. These tiny forces are hardly detectable even with very sensitive 
methods such as atomic force microscopy but nonetheless they are still strong enough to affect 
the Brownian motion of particles of micrometer size (i.e., colloids) which can be used as detectors. 
This suggests the study of a geometrical arrangement in which the critical fluctuations of the 
medium are restricted in the space between a spherical colloid and a planar surface rather than 
between two parallel surfaces, as considered so far. On the basis of the discussion of subsec. \T2\ 
one actually expects a Casimir force F to act on the sphere. If the minimal distance z between 
the surface of the sphere and the planar substrate is much smaller than the radius R of the 
sphere itself, the spherical surface can be approximated by a set of circular rings parallel to the 
opposing planar substrate (Derjaguin approximation) and therefore the force in this geometrical 
setting can be related to the one for parallel plates. Within this approximation, the z-dependent 
critical Casimir force F at temperature T acting on the sphere is given by (see ref. [2^ for details) 



where the universal scaling function -d^o can be calculated in terms of a suitable integral of the 
scaling function of a film [25] and it retains the qualitative features of i^n. The theoretical 
prediction for the scaling form of the Casimir potential ^c{z) = dz'F{z') = k^T {R / z)'E{z / S^) 
immediately follows from the integration of eq. ([H]) . 

In order to measure tiny forces acting on a single colloid, one can take advantage of the Total 
Internal Reflection Microscopy (TIRM) [23] , a very sensitive technique that is widely employed 
in the study of soft-matter systems |23j . With TIRM one monitors the Brownian motion of 
the colloid floating in a fluid and from the associated statistics one infers the potential of the 
forces acting on the particle. The scheme of the experimental setup is depicted in fig. [5]^a): 
an incoming visible laser beam is totally reflected at the substrate (glass)-fluid interface, so 
that an evanescent field of intensity Iew{xA_) penetrates the liquid with a typical exponential 
decay as a function of the distance x± from the interface, see fig. [SU^b). The colloidal particle, 
of suitable refractive index and at a surface-to-surface distance z from the substrate, scatters 
light out of this evanescent field. The intensity I^c of the scattered light is measured by a 
photomultiplier (PM) and turns out to be proportional to the intensity Iev{z) of the evanescent 
field. From the time dependence of the scattered intensity Isc{t) oc /ev(-z(i))j represented in 
fig- [Mc)) one can therefore calculate the time dependence of the particle- wall separation, i.e., 
z{t). This dependence, depicted in fig. [5]|^d), refiects the Brownian motion of the particle under 
the effects of thermal fiuctuations and possible additional forces. The probability distribution 
P{z)., as extracted from the histogram of {•z(i)}o<t<tsamp for sufficiently large sampling time 
isamp (fig- [M6))i is proportional to the Boltzmann factor exp{— /3$(2:)} (with [5 = {k^T)'"^), 
where ^{z) is the potential of the total force acting on the particle. Accordingly, $(z) can be 
calculated up to an irrelevant constant on the basis of the P{z) inferred from the experimental 
data. 




(8) 




Figure 5. Determination of the potential ^(z) of the forces acting on a single colloidal particle 
by means of the Total Internal Reflection Microscopy. 



In order to provide a direct evidence of critical Casimir forces one can study a single colloidal 
particle immersed in a binary liquid mixture of water (W) and lutidine (L), close to a glass 
substrate [25]. In ref. [25j TIRM has been employed for the determination of the potential $(z) 
of the forces acting on the particle and for studying how $ changes when the binary mixture is 
driven towards the critical point. The W-L mixture has a bulk phase diagram such as the one 
presented in fig.Hl with a lower critical point at Tc = 34°C and lutidine mass fraction cl,c — 0.29. 
The local excess concentration of lutidine (5cl (x) is a suitable order parameter for this demixing 
transition. In what follows we illustrate the results of the experiments in the case of a mixture 
at the critical concentration cl,c (the case cl 7^ cl,c is discussed in refs. [Ml |25] ). 

At a temperature T far below the critical value Tc, fluctuations of the order parameter 
are irrelevant and the potential ^(z) of the forces acting on the particle is basically the 
sum of two contributions: (a) a screened electrostatic repulsion (<I>ei(-z) oc e"**^) between the 
colloid and the glass substrate, relevant in our experimental conditions for z < 0.1 /im and 
(b) a linearly increasing potential due to the combined effect of buoyancy, gravity and optical 
pressures generated by the optical tweezer employed in the actual experimental setupH. This 
latter contribution has been subtracted from the curves presented in fig. El which show the 
dependence of the measured potential $(z) on the distance from the critical point. Panel (a) 

and (b) in fig. prefer to (colloid, substrate) = ( ) and (H — ) BCs, respectively. As expected, 

few hundreds mK far from the experimentally estimated critical temperature, the potential ^{z) 
consists only of the electrostatic repulsion, which is negligible on this scale for z > 0.12 ^m in 
(a) and 2; > 0.1 ^m in (b). Upon increasing the temperature by few tens mK significant changes 
occur in the measured potentials as a strong attractive and repulsive force develops, respectively, 
in (a) and (b). In particular, in both cases one observes a change in the total potential of the 
order of 10 k^T, which is really significant on the scale of ksT which characterizes the physics 
of colloids. The maximum value of the attractive force acting in case (a) is estimated at 600 
fN. The strong temperature dependence of the measured potentials is a clear indication of the 
involvement of critical Casimir forces. This evidence is supported by the comparison with the 



^ A third possible contribution, due to van-der-Waals forces turned out to be negligible in the experiment reported 
in ref. [^[^ . 




Figure 6. Temperature dependence of the potential ^{z) of the forces acting on single 
polystyrene colloidal particles at a distance z of closest approach from a glass substrate. 
The particles and the substrate are immersed in a water-lutidine liquid mixture with critical 
composition and temperature T smaller than the critical one Tc. Data points refer to the 
experimental results (where gravity and buoyancy have been subtracted, see the main text) 
obtained for a substrate treated with NaOH in order to make it hydrophilic [(— ) BC for the 
order parameter] and (a) a colloid of diameter 2R = 2.4 /Ltm absorbing preferentially water (— ) 
(b) a colloid of diameter 2R = 3.7 //m with preferential adsorption for lutidine (+). The solid 
lines in (a) and (b) are the corresponding theoretical predictions for the critical Casimir potential 
with (colloid, substrate) = ( ) and (H — ) BCs, respectively. (From ref. p5].) 



corresponding theoretical predictions for <I>c'(z), shown as solid lines in fig. O In order to focus 
only on the contribution to ^{z) due to the critical Casimir force, this comparison is presented 
only within the range of distances where the electrostatic contribution is negligible. The limited 
accuracy with which the critical temperature Tc has been experimentally determined in ref. [25j 
does not allow a reliable calculation of the correlation length ^ via ^ = ^o(l ~ T/Tc)~^'^^ and 
on the basis of the value ^q*^^^^ of measured for the W-L mixture in the bulk by independent 
experiments [26] . Instead, for each single temperature T, the correlation length (,(T) can be 
determined in order to optimize the agreement between theory and experimental data. Then, 
separately for the two experiments, the value ^g^*^ of and Tc^*^ of Tc are determined in such a 
way to yield the best fit of the determined (,(T) with the expected algebraic behaviour. In both 
cases the fact that, within errorbars, ^q^*'* = ^o''^^^ is a check of the significant agreement between 
theory and experiments. (See refs. |241I25| for a more detailed discussion.) Interestingly enough, 
via a suitable chemical surface treatment, it is possible to change quite easily the substrate from 
hydrophilic (— ) to hydrophobic (+) and therefore switch the Casimir force acting on the colloids 
in fig. [6l||a) and (b), respectively, from attractive to repulsive and vice-versa (see refs. [MIES]). 

Summing up, by measuring their effect on the Brownian motion of a colloidal particle, direct 
evidence of both attractive and repulsive critical Casimir forces has been provided. Based on 
the general arguments presented in subsec. 12.21 we expect such forces to act also between two 
or more colloids immersed in a near-critical mixture. Due to the strong non-additivity of these 
fluctuation-induced forces we expect interesting many-body effects. 



4. Conclusions, perspectives and applications 

4-1. Conclusions 

In the previous sections we presented an overview of the Casimir effect due to critical fluctuations, 
highlighting its relation with the analogous effect in QED. We focussed on the theoretical 
predictions which are relevant for the interpretation of the available and possibly forthcoming 
experiments, especially involving wetting films and colloidal particles. The comparison between 
the available theoretical and experimental results turns out to be very good. We did not attempt, 
however, to review of all the currently available relevant theoretical results and approaches, a 
task which is well beyond the scope of the present contribution. In what follows we mention 
only some of them. 

The discussion in subsecESl clearly shows that Casimir-like effects are expected whenever the 
spectrum of fluctuations of a certain nature is changed by the presence of confining bodies, which 
consequently experience fluctuation-induced forces. The range of such forces is typically set by 
the correlation length ^ of the relevant fluctuations and a universal behaviour has generally to 
be expected whenever ^, L » ^micr where L is the typical distance between the boundaries and 
imicr a microscopic, system-dependent scale. 

4-2. Perspectives 

In view of the quite general conditions under which such Casimir forces arise, several examples 
can be found in statistical physics. They occur not only at the critical points previously 
discussed, or in systems belonging to other bulk and surface universality classes (e.g., liquid 
crystals or complex fluids [71 EZ]), but also close to Uncritical points. This instance has been 
studied both experimentally and theoretically [28l|T7j. Due to the fact that the upper critical 
spatial dimensionaliy for a tricritcal point is three, a mean-field approximation of suitable models 
captures rather well the features observed in experiments. Casimir-like forces result also from 
(long-ranged) fluctuations occurring in non-equilibrium steady states of different nature, from 
chemical reactions to granular matter, see, e.g., ref. [221, or from quantum fluctuations of the 
order parameter of a quantum phase transition [llj. However, at present it is not clear how 
these predictions can be experimentally tested. On the other hand, there are experimentally 
relevant systems, such as Bose gases close to the Bose-Einstein condensation, which have not yet 
been studied with the aim of detecting critical Casimir forces, even though some quantitative 
theoretical predictions (see fig. [3l^b)) are available beyond the case of ideal gases |30] . Thermal 
capillary fluctuations of the interface between two different liquids give rise to a Casimir force 
acting on two or more collodal particles trapped at that interface [31]. In this case, the presence 
of the colloids alters the spectrum of these surface fluctuations, with the additional feature that 
the boundaries themselves at which the boundary conditions are imposed, represented by the 
contact line of the interface at the surface of the colloids, fluctuate (see ref. [31j for details). The 
resulting forces might be of relevance for understanding the topical problem of the interaction 
among colloids trapped at interfaces. 

So far we have only considered homogeneous substrates and particles. However, in the 
presence of chemical or topographical modulations, such as the one depicted in fig. [TJ novel 
phenomena are expected, including the emergence of lateral critical Casimir forces and torques 
Tcas! still characterized by a universal scaling behaviour. Some of the features of these forces 
have been studied in refs. [32 j and might form the basis for the quantitative understanding of 
recent experimental investigation of the behaviour of a single colloidal particle immersed in a 
binary mixture and exposed to a patterned substrate [33] . 

In addition to the equilibrium static behaviour of the critical Casimir effect, its equilibrium 
and non-equilibrium dynamics poses both a theoretical and experimental challenge. Indeed, it 
is not even obvious how to define theoretically such forces when their definition based on the 
decomposition of the effective free energy (see eq. ([5]) ) is no longer viable due to the lack of a 




(a) (b) (c) (d) 

Figure 7. Possible chemically [(a), (d)] and topographically [(b)] patterned substrates, in which 
a lateral critical Casimir force -Fy arises in addition to the normal force F± previously discussed. 
A critical Casimir torque rcas is expected to act on non-spherical particles [(c)]. 

proper notion of free energy for dynamical phenomena. A rich dynamical behaviour is actually 
emerging already in quite simple theoretical relaxational models |34j . In order to make contact 
with the dynamical information which might come from TIRM experiments, it remains to be 
seen how the time-dependence of the critical Casimir force affects the Brownian motion of the 
colloidal particle in conjunction with hydrodynamic and adsorption phenomena. 

4-3. Applications? 

In contrast with the interactions typically acting among colloids (e.g., elctrostatic) , the critical 
Casimir forces show a striking temperature dependence, as we have reported in subsec. [3?2l This 
fact can possibly be exploited in order to control via minute temperature changes the phase 
behaviour and aggregation phenomena in systems with dispersed colloids. In addition, not only 
the range of interaction can be controlled but also the sign of the resulting force. Typically 
this can be achieved by surface treatments and does not require (as it does for the quantum 
mechanical Casimir effect) substantial changes or tuning of the properties of the bulk materials 
which constitutes the immersed objects. This property might be exploited in order to neutralize 
the attractive quantum mechanical Casimir force responsible for the stiction which brings micro- 
electromecanical systems to a standstill. If these machines would work not in a vacuum but in 
a liquid mixture close to the critical point, the stiction could be prevented by tuning the critical 
Casimir force to be repulsive via a suitable coating of the various machine parts. In principle, 
by using optically removable or controllable coatings, one could very conviniently control the 
functioning of the microdevice without acting directly on it. 
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